Isolated Vibrational Wavepackets in D^: Defining Superposition Conditions and 

Wavepacket Distinguishability 
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Tunnel ionization of room-temperature D2 in an ultrashort (12 femtosecond) near infra-red (800 
nm) pump laser pulse excites a vibrational wavepacket in the ions; a rotational wavepacket is 
also excited in residual D2 molecules. Both wavepacket types are collapsed a variable time later 
by an ultrashort probe pulse. We isolate the vibrational wavepacket and quantify its evolution 
dynamics through theoretical comparison. Requirements for quantum computation (initial coher- 
ence and quantum state retrieval) are studied using this well-defined (small number of initial states 
at room temperature, initial wavepacket spatially localized) single-electron molecular prototype by 
temporally stretching the pump and probe pulses. 

PACS numbers: 42.50.Hz, 82.53.Hn, 03.67.Lx 
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As with a broad spectrum of quantum systems, if a 
molecule laser-induced into a coherent superposition of 
vibrational states is free to evolve, a time-dependent 
wavepacket is created. The wavepacket spreads or de- 
phases significantly after the initial superposition, the 
result of the constraining potential surface anharmonic- 
ity. However, when the discrete eigenstates naturally 
rephase or revive, the original state of the system is al- 
most fully recovered, and the initial quasi-classical peri- 
odicity returns. This phenomenon of wavepacket revival 
[H has major implications in experimentally realisable 
time-evolving quantum systems. Through the precise op- 
tical preparation of a tailored wavepacket in a molecule, 
the wavepacket motion is imprinted on the dissociative 
states, observable with a probability governed by the am- 
plitudes that contribute to the superposition. 

In this Letter, we report the isolation and quantifica- 
tion of a vibrational wavepacket in a room-temperature 
ensemble of ions. We vary the initial coherence of 
the superposed states and distinguish structure in the 
wavepacket revival when the measurement condition is 
changed. An application of these results is as a single- 
internal degree-of-freedom prototype vibrational qubit 
(vibrational states, v, support the qubit), the building 
block of a molecular quantum computer (MQC). Quan- 
tum computation operations are not performed, rather 
we test the preparation and observation of a superposi- 
tion. While MQC has been demonstrated in Li2 the 
system is considered as it is readily handled theoret- 
ically (numerical solution of the rotation-free ID TDSE 
under the two-state approximation) Q and generated ex- 
perimentally (D2 molecules at room temperature are pre- 
dominantly in u = 0, no state-selective cooling is required; 
short laser pulse defines a spatially localized initial state) 



Before implementing quantum computation (QC), 
five general 'physical implementation' requirements have 
been identified Q: (1) the ability to initialize the qubit 
in a fiducial state, (2) the coherence time of the qubit is 
much longer than the operation times, (3) a universal set 
of logic operators are available, (4) a qubit-specific mea- 
surement capability exists and (5) the quantum system 
is well characterized and readily scalable. We quantify 
the fulfilment of requirements (1), (4) and (5) with a dis- 
cussion of requirement (2). 

Clearly requirement (3) lies at the heart of QC, which 
in future experiments will be attempted by optically con- 
trolling the phase 0] or amplitude Q of each vibrational 
state, referred to as phase or amplitude manipulation. 
The additional degrees of freedom in a larger molecule 
could be employed i.e. compound vibrational or 
bending modes, making scalable (multi-dimensional) QC 
feasible. 

Tunnel ionization (D2 {v' = 0) -I- pump — > {v = 
0, 1, 2...) -I- e~) when the Keldysh parameter 7 < 1 
creates a vibrational wavepacket on the IsCTg potential 
surface of [l, 0, HI, and is near- instantaneous so is 
treated as a 'vertical' transition. The configuration of 
the potential energy surfaces is such that ionization only 
transfers population around the inner turning point of 
the wavepacket motion 0, 0, S] ■ In a phase manipulation 
MQC, such a transition fulfils requirement (1), as all v- 
states are initiated with essentially zero phase difference. 
In an amplitude manipulation MQC, requirement (1) is 
similarly fulfilled as all accessible w-states have a known 
(Franck-Condon) amplitude distribution However, 
in practice an ensemble of molecules may be used to gen- 
erate meaningful signal, where it becomes important that 
the pulse length is sufficiently short to retain coherence 
in the ensemble. 
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FIG. 1: (Color online) Comparison between experiment and 
theory showing vibrational and rotational revivals, (a) 
PD (top) and DI (bottom) yield as a function of wavepacket 
evolution time, At, and ion kinetic energy, recorded at room 
temperature, (b) Integrated PD yield (thick black) contains 
vibrational and rotational revivals (top) . The theoretical pre- 
diction (green / thick grey curve) shows remarkable agree- 
ment with the experimental results when rotational dynam- 
ics are excluded (thin black curve), (c) The calculated time- 
dependent probability distribution, (d) The calculated impul- 
sive alignment of D2 induced by the pump pulse in a room- 
temperature distribution of rotational levels, (e) The pre- 
dicted Franck-Condon (FC) distribution, (f) Fourier Trans- 
form of the experimental PD yield. The frequency compo- 
nents present in the vibrational wavepacket are beats between 
vibrational states v and v ± 1 . The Fourier amplitudes thus 
return a measure of vibrational population. 



Observing the vibrational wavepacket evolution on the 
IsCTg ground potential surface is possible by 'collapsing' 
the wavepacket to an observable state with a second laser 
pulse. Photodissociation (PD) -I- probe — > D"*" -I- D 
breaks the internuclear bond via photon coupling the 

ISCTo 



sec- 



/g ground and 2pcr„ first excited states 
ond route to observing wavepacket evolution is through 
double ionization (DI) DJ D+ + D+ + e^, requiring 
a higher laser intensity as a second ionization event must 
occur I4I . PD and DI are distinguishable by the kinetic 
energy released during bond fracture. Both PD and DI 
are enhanced when the wavepacket is at the outer turning 
point ~ 3.5 a.u., [3, and the internuclear potential 
energy surfaces naturally act as a quantum 'shutter' or 
'barrier' depending on whether they are laser-coupled or 
not. 



The experimental set-up is an extension of one de- 
scribed elsewhere 0, i, [l|: intense (4 x lO" Wcm'^) 
ultrashort (12 femtosecond, 1 fs 10^^^ s) near infra- 
red (800 nm) laser pulses are used to create (pump) and 
observe (probe) the vibrational wavepacket. A co-linear 
low-dispersion Mach-Zehnder interferometer capable of 
supporting the pulse bandwidth (> 100 nm) generates 
the time-separated pump and probe, allowing a delay 
< At < 160 picoseconds between the pump and probe 
with a repeatable resolution of ~ 300 attosecond (1 as = 
10~^^ s). An alternative method for generating two or 
three time-delayed pulses has been demonstrated [13] by 
transmitting a ~ 10 fs pulse through a disc-double annu- 
lus optic allowing two variable delays of hundreds of fs 
if the two annuli are radially rotated. After focussing 
with a spherical mirror, the interaction volumes gen- 
erated by our Mach-Zehnder configuration overlap per- 
fectly, whereas in the disc-annulus configuration, a time- 
varying spatial overlap is produced. Furthermore, the 
dispersion of the Mach-Zehnder configuration does not 
change with delay. D+ fragments from PD and DI, Fig. 
1(a), are detected in an ion time-of- flight mass spectrom- 
eter containing a 250/xm aperture, thus signal can only 
result from PD or DI when the molecular axis is parallel 
to the probe polarization and detector axis. The aper- 
ture also defines a quasi-isointensity source in the focal 
volume. 

In addition to creating a D^ vibrational wavepacket, 
few-cycle pulses impulsively align deuterium [l3|, Il4| . 



thus the revival of a rotational wavepacket is expected. 
We isolate the vibrational and rotational motions, which 
is only possible for a limited range of pump pulse dura- 
tions as discussed later. As shown in Fig. 1(a), in the 
region 450 < At < 650 fs, both PD and DI yields exhibit 
a ~ 24 fs revival modulation with maximum contrast at 
550 fs, the result of vibrational revival, consistent with 
recent findings [3, 0, S]- However, between 250 < At < 
350 fs a rotational revival is also observed, comparable to 
recent observations (isl . [l3 | which will also underlie the 
vibrational revival around 550 fs. Both the rotational 
and vibrational revivals are evident in the integrated PD 
yield as a function of At, Fig. 1(b) (top). 

A major finding of the present work concerns the sep- 
aration of these two wavepackets. The isolation of the 
vibrational wavepacket requires confirmation that two 
distinct populations in D2 and D^ exist followin g th e 
pump. In recent work by the authors [3] (see also [3]), 
we quantify the rotational wavepacket in D2 by Fourier 
transformation of the high-energy time-dependent DI sig- 
nal recorded under identical experimental conditions. We 
observe peaks (with an accuracy of 5%) in the frequency 
spectrum at 5.5 (0), 9.2 (1), 12.8 (2), 16.4 (3), 20.1 (4) 
and 23.7 (5) THz, where parenthesis indicate the J : J -I- 
2 rotational state beat in the D2 molecule. There is no 
evidence for a rotational wavepacket in D J , which would 
contribute peaks in the frequency spectrum at 2.7 (0), 
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FIG. 2: (Color online) Controlling the initial superposition 
in an isolated vibrational wavepacket. (a) Experimental PD 
yield as the pump duration is increased (thin black curve). 
The rotational wavepacket has been removed from the dataset 
as discussed. The shape of the revival indicates the superpo- 
sition degradation as the pump is stretched, illustrated by 
comparison with a theoretical prediction (thick magenta / 
grey curves), (b) The predicted vibrational wavepacket after 
averaging over the changing Gaussian pump duration. 



4.5 (1), 6.3 (2), 8.1 (3), 9.9 (4), 11.7 (5), 13.5 (6), 15.3 
(7) and 17.1 (8) THz. This is confirmed in the present 
work. Furthermore, the vibrational wavepacket can only 
exist in the molecular ion generated in the pump, 
as the pump duration is greater than the quasi-classical 
period (11 fs) of D2 as recently investigated 15]. 

Thus we find that the pump either creates a vi- 
brational wavepacket in (rotating in a non-time- 
correlated manner with respect to the pump with a 
room-temperature distribution of J-states) or a rota- 
tional wavepacket in D2 (vibrating incoherently in an im- 
pulsively modified distribution of J-states). Solving the 
TDSE for the vibrational (Fig. 1(c)), and rotational 
[1^ (Fig. 1(d)) degrees of freedom, we quantify these 
motions. We take the integrated PD yield as a function 
of At (Fig. 1(b) top), and subtract the probability of 
finding the rotational wavepacket within the detector ac- 
ceptance. Fig 1(d). The integrated PD yield (arb. units) 
is scaled relative to the rotational wavepacket probability 
until the signature of the rotational wavepacket (250 < 
At < 350 fs) is totally suppressed. The resulting PD yield 
contains the pure vibrational wavepacket, as shown in 
Fig. 1(b) (bottom). The PD yield is modelled by numer- 
ically propagating the initially localized state (FC princi- 
ple) then calculating the time dependent coupling of the 
IsCTg and 2pcr„ states The vibrational wavepacket in 
the molecular ion is clearly well-described by this rela- 
tively simple treatment, thus requirement (5) is partially 
fulfilled. If evidence for strong rotational- vibrational cou- 
pling had been found, this would not be the case. 

Fourier transformation of the raw time-dependant PD 
yield. Fig 1(b), recovers the beat frequencies of the vi- 
brational levels with the associated amplitudes that form 



the vibrational wavepacket, therefore we have a method 
for recovering the population of the vibrational states. 
As is apparent from Fig. 1(e) and (f), the measured 
distribution of vibrational beat frequencies agrees with 
the expected FC distribution of states for w > 4. For 
states V < 4, negligible beating is observed. This is a 
result of the probe pulse being of insufficient intensity to 
open the three-photon crossing 12], thus the vibrational 
wavepacket 'below' v = 4 remains bound. Our 'quantum 
shutter' is therefore selective, however the principle is 
sound, since given a probe intensity > 6 x 10^^ Wcm~^ 
all V states will be accessed. It should be noted that the 
states < ?; < 4 must be populated, otherwise the quasi- 
classical period of Dj would not be observed: lack of 
population in these states would increase the vibrational 
revival periodicity. In a phase manipulation MQC, the 
relative phase can be temporally resolved for occupied 
Instates; while in an amplitude manipulation MQC, the 
relative amplitudes of the constituent instates can also be 
recovered by Fourier transform analysis, Fig. 1(f). Re- 
turning to QC requirement (4), the ability to determine 
the final state in our measurement is met. Destruction of 
the prototype qubits is necessary to read the final state, 
however it is not necessary to destroy the wavepacket to 
manipulate. By utilizing nondestructive pulses, it should 
be possible to modify the wavepacket while maintaining 
the superposition. 

To extend our understanding of fulfilling requirements 
(1) and (4), we quantify variations in the initial vibra- 
tional wavepacket coherence. We also study wavepacket 
distinguishability, the faculty to temporally discern sub- 
sequent oscillations of the vibrational wavepacket, partic- 
ularly around revival. Both are macroscopically achieved 
through independently varying the pump and probe du- 
rations by changing the dispersion (hence chirp) in either 
arm of the interferometer. 

Fig. 2(a) shows the integrated PD yield in the vicin- 
ity of the vibrational revival as the pump duration, rp„, 
is increased from 14 to 28 fs while the probe duration, 
Tpr, is fixed at 13 fs. The vibrational wavepacket in DJ 
was isolated from the rotational signature in D2 as de- 
scribed earlier. The contrast of the vibrational revival 
structure decreases as the pump duration increases, the 
consequence of loss of coherence in the initial super- 
posed state. The experimental revival structure (thin 
black curve. Fig. 2(a)) is compared to the theoretical 
wavepacket (Fig. 1(c)) averaged over the pump pulse. 
Fig. 2(b), and integrated 3 < R < 10 a.u., then over- 
laid (thick magenta / grey curves) on the experimental 
data in Fig. 2(a). The resultant magenta curves show 
an excellent agreement with the decohering isolated vi- 
brational wavepacket, clearly demonstrating the role of 
coherence in the initial pump pulse. 

In Fig. 2(a), it is apparent that even for Tp„ — 28 fs, 
some vibrational coherence remains when Tp„ is longer 
than the revival oscillatory period, 24fs. This residual 
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FIG. 3: (Color online) Distinguishing the reviving vibrational 
wavepacket as the measurement condition is changed, (a) The 
sudden loss of revival fidelity when the probe duration exceeds 
40 fs, the consequence of the temporal structure becoming in- 
distinguishable to the probe. From top to bottom, the probe 
durations are indicated sequentially. Comparing a (b) 15 fs, 
(c) 30 fs and (d) 45 fs probe to the revival structure, the 
facility to decern subsequent oscillations (i.e. wavepacket dis- 
tinguishability) of the vibrational revival is compromised. 



structure is due to the highest vibrational states, v > 9 
which have the longest periods and an absolute popula- 
tion of less than 4%. Consequently, care must be taken in 
defining the initial state of a vibrational qubit as states 
with an apparently negligible population clearly influence 
wavepacket evolution. 

By increasing Tpr from 14 to 60 fs while keeping 
the pump fixed at Tp„ = 14 fs (thus the vibrational 
wavepacket is identical to Fig. 2(b), bottom, through- 
out), the amplitude of the revival changes suddenly as 
Tpr exceeds 40 fs. Since a 15 fs probe. Fig. 3(b), is short 
compared to the revival, oscillations are observed with 
high fidelity. At 30 fs. Fig. 3(c), a comparable fidelity 
is also achieved, as subsequent revivals are still distin- 
guished. However, a 45 fs probe. Fig. 3(d), cannot dis- 
tinguish any two revival maxima, hence the oscillatory 
structure is suddenly lost when Tpr approaches double 
the quasi-classical period of . We can therefore define 
a wavepacket distingushability condition, which is less 
strict than for creation, depending on the measurement 
regime. If a time-dependent measurement is performed 
at revival (as required in the phase manipulation), Tpr 
must be shorter than double the quasi-classical period to 
ensure total wavepacket collapse. Alternatively, to guar- 
antee that all the wavepacket is measured in an energy- 
resolved scheme (as required in the amplitude manipu- 
lation), Tpr of at least double the period of the highest 
populated state is required. Furthermore, the intensity 
requirements discussed earlier must also be adhered to. 

The pump-variation observations indicate a promis- 
ing future for controlling the initial superposition con- 



ditions. Following the optical preparation of an iso- 
lated vibrational wavepacket in a single molecule, a ge- 
netic algorithm will determine the intermediate time- and 
frequency-structured laser pulse [l6| required to achieve 
a particular target state (either by phase or amplitude 
manipulation). This intermediate 'operator' pulse must 
be low intensity so as not to fracture the internuclear 
bonds, rather to guide wavepacket evolution. The high 
intensity ultrashort probe pulse will then fracture the in- 
ternuclear bonds, allowing the result of the operation to 
be read. 

To perform logic operations and apply successful error 
correction schemes requires a coherence time ~ IC* times 
the interaction time, t/, requirement (2). Here, ioniza- 
tion defines t/ ~ 1.3 fs thus if the wavepacket survives 
to ~ 10 ps with high fidelity, QC is possible. Our mea- 
surements have exceeded 2 ps, therefore this constraint 
does not appear technologically limiting. The MQC will 
be scaled by increasing the number of vibrational modes 
supporting the qubit, accomplished in a larger molecule. 
Communication between states will then be controlled 
optically by distorting the electronic potential surface 
through the Stark-shift. 

In a broader application, a pump of a shorter duration 
than the vibrational motion establishes a time-dependent 
variation of the refractive index in a Raman-active molec- 
ular medium, as demonstrated in SFg 17]. By employing 
the vibrational mode of a hydrogenic molecular ion, it 
may be possible to spectrally broaden and self-compress, 
potentially to the single-cycle level. 
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